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ABSTRACT 

Cosmological simulations show that, at the present time, a substantial fraction of the gas in the intergalactic 
medium (IGM) has been shock-heated to T > 10 5 K. Here we develop an analytic model to describ e the frac- 
tion of shocked, moderately overdense gas in the IGM. The model is an extension of the iPress & SchechteJ 
i 1974) description for the mass function of halos: we assume that large-scale structure shocks occur at a fixed 
overdensity during nonlinear collapse. This in turn allows us to compute the fraction of gas at a given redshift 
that has been shock-heated to a specified temperature. We show that, if strong shocks occur at turnaround, 
our model provides a reasonable description of the temperature distribution seen in cosmological simulations 
at z ~ 0, although it does overestimate the importance of weak shocks. We then apply our model to shocks at 
high redshifts. We show that, before reionization, the thermal energy of the IGM is dominated by large-scale 
structure shocks (rather than virialized objects). These shocks can have a variety of effects, including strip- 
ping < 10% of the gas from dark matter minihalos, accelerating cosmic rays, and creating a diffuse radiation 
background from inverse Compton and cooling radiation. This radiation background develops before the first 
stars form and could have measurable effects on molecular hydrogen formation and the spin temperature of the 
21 cm transition of neutral hydrogen. Finally, we show that shock-heating will also be directly detectable by 
redshifted 21 cm measurements of the neutral IGM in the young universe. 

Subject headings: cosmology: theory - intergalactic medium - large-scale structure of the universe 



1. INTRODUCTION 

One of the major achievements of cosmology over the past 
two decades has been the development of a coherent and 
testable paradigm for structure formation. In the modern pic- 
ture, small density inhomogeneities imprinted on the universe 
during inflation grow continuously through gravitational in- 
stability. These density fluctuations still have small amplitude 
at the time of cosmological recombination (z ~ 1000), when 
the cosmic microwave background (CMB) is last scattered. 
At this epoch, structure formation is straightforward to study 
analytically because the fluctuations are still linear. Observa- 
tions of the CMB have shown convincing evidence that our 
theories of structure formatio n accurately descri be the uni- 
verse at this early phase (e.g., Spergel et al. 2003). Eventu- 
ally, the perturbation amplitude grows sufficiently large for 
nonlinear effects to become important: at this point the first 
gravitationally bound structures begin to form. The struc- 
ture formation problem then passes beyond the scope of ana- 
lytic approache s (except in the case of perfectly spherical col- 
lapse; Gunn & Gott 1972). The most successful approxima- 
tion for studying the properties of gravitational instability in 
the n onlinear regime is the so-called Zel'dovich approxima- 
tion llZeF Do vichl l 1 9701) . which uses Lagrangian (rather than 
Eulerian) perturbation theory. The theory shows that collapse 
will occur first along a single axis, creating a "pancake" or 
sheet, then along a second axis, creating a filament, before 
finally converging towards a single point. 

Fully non-linear treatments of structure formation had to 
await the development of numerical simulations. These sim- 
ulations confirm the behavior predicted by the Zel'dovich 
approximation: at low redshifts, the matter in the universe 
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organizes itself into a complex network of sheets and fila- 
ments (with galaxies and galaxy clusters at their intersec- 
tions) known as the "cosmic web." This paradigm has had 
admirable success reproducing the large -scale distribution of 
galaxies in redshift surveys (e.g.. Ide L apparent et al. 1986) 
and in describi ng the absorption features of the Lyq for- 
est (e .g., ICen et al Jfl99l iHernquist et alJll996tlTheuns et alJ 
1998). 

While the dark matter accretes smoothly onto the cosmic 
web, the baryons do not: their infall velocities onto these 
structures often exceed the local sound speed and a com- 
plex network of shocks forms. Probably most important 
are the well-known virial shocks that form when halos col- 
lapse. However, cosmological simulations also contain many 
shocks outside of collapsed halos; th ese fill a fair fraction 
(~ 10%) of the volume of the univer se llCen & Ostrike rTl999t 
iKeshet et al1l2003t iRvu et alJ 120031) and play an important 
role in shaping the properties of baryons in the local uni- 
verse. Simulations predict that > 20% of all baryons are in 
the "w arm-h ot intergalactic medium" (WHIM) with T ~ 10 5 - 
10 7 K (Dave et all2001l) . most of which have densities below 
those characteristic of virialized objects. 

Unfortunately, there is little analytic understanding of the 
distributio n of shocks in the intergalactic medium (IGM). 
iCen & OstrikeJ dl999l) argued that the characteristic temper- 
ature of the gas can be computed by simply associating the 
postshock temperature with a fluctuation collapsing over the 
age of the universe that has just become nonlinear at the 
present epoch. This provides a rough guide to the expected 
strength of the shocks but does not predict how much gas 
has been shocked or give a detailed temperature distribution. 
iNath & Trenthamj {T997) expanded this treatment using the 
Zel'dovich approximation to give rough estimates of the frac- 
tion of gas in the WHIM phase, but they were not able to 
predict the detailed properties of shocks at any redshift. 

There has also been little consideration of the role that sim- 
ilar large-scale structure shocks play at higher redshifts. One 
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reason that they are significant in the local univers e is that the 
char acteristic shock temperature (r s h ~ 10 7 K; see Dave et al. 
1200 ll) is much larger than the mean temperature of the pho- 
toionized IGM (7igm ~ 10 4 K). As redshift increases, de- 
cre ases rapidly; at z > 2, t he characteristic temperature given 
by ICen & Ostrikerl (1999) is not much larger than 7igm, so 
we expect strong shocks to be much less common. How- 
ever, before reionization, the typical temperature of the IGM 
is much smaller and shocks can again become important (de- 
spite the much smaller-scale perturbations that collapse). Be- 
cause shocks occur in overdense regions, where galaxies also 
preferentially occur, they obviously help to determine the en- 
vironments in which most galaxies reside. The large-scale 
environment has important consequences for feedback from 
(and on) collapsed objects, but this point has not been consid- 
ered in the past. 

In this paper, we consider a new analytic model for the tem- 
perature distribution of shocked, moderately overdense gas in 
the IGM. We are guided by the fact that the distribution of 
collapsed objects is well-described, to a surprisingly high pre- 
cision, by the Press-Schechter formalism ( Press & Schechter 
[1971 . The key ingredients to this model are (1) Gaussian 
statistics for the density field and (2) spherical collapse of 
perturbations in the nonlinear regime. Motivated by the suc- 
cess of this formalism, we use the same assumptions to model 
shocked gas in the IGM. We show that our simple model 
provides a reasonable fit to existing (low-redshift) simulation 
data if only strong shocks are considered. In |2] we present 
our formalism for the large-scale structure shocks and use it to 
describe the basic properties of shocks at both low and high 
redshifts. We then move on to consider the role that high- 
redshift shocks might play. In ^3] we describe three impor- 
tant effects that the shocks can have: accelerating particles to 
high energies, creating a diffuse background radiation field 
(principally through cooling radiation), and stripping mini- 
halo gas. In J|we consider the observable consequences of 
high-redshift shocks. Finally, we summarize our main con- 
clusions in 50 

Throughout the discussion we assume a flat, A-dominated 
cosmology, with density parameters fi m = 0.3, f^A = 0.7, and 
£lb = 0.047 in matter, cosmological constant, and baryons, re- 
spectively. In the numerical calculations, we assume h = 0.7, 
where the Hubble constant is Ho = 100/z km s -1 Mpc" 1 . We 
take a scale-invariant primordial power spectrum (« = 1) nor- 
malized to an amplitude of er^ = 0.9 on 8/T 1 Mpc spheres. 
These values are co nsistent with the most recent measure- 
ments dSnergel et al 120031). 

2. A MODEL OF SHOCKS FROM LARGE-SCALE STRUCTURE 

Following lCen & Ostrikerl (1 19991) and lKeshet et alJ J2003I) . 
we associate large-scale shocks with the collapse of nonlinear 
objects. The shock velocity v s is then approximately equal to 
the length scale of the perturbation divided by the age of the 
universe, or 

v s ^H(z)R pi (1) 

where R p is the physical radius that the region would have had 
if it expanded uniformly with the Hubble flow (i.e., 4nR p /3 = 
M s /p z , where M s is the mass of the shocked region and p z is 
the mean density of the universe at redshift z). The postshock 
temperature T s h is then 
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where fjan p is the mean particle mass, 7 is the adiabatic in- 
dex of the gas, and r\ is the square of the Mach number of the 
shock in the downstream medium. Given a shock velocity and 
the sound speed of the IGM gas, r\ is fixed in terms of 7 by the 
jump conditions; however, the definition of the shock velocity 
in equation Q is only approximate, and we let r\ encapsulate 
uncertainties in both the collapse time and the relevant length 
scale. In most of the following discussion we will assume 
strong shocks, for which r\ = 5/16 (if 7 = 5/3). Note that the 
shock temperature depends on both th e mass of t he perturba- 
tion an d on redshift. ICen & Ostrikerl (|l999) and lDave et alJ 
(2001) argued that choosing R p to be the length scale associ- 
ated with the characteristic nonlinear mass at a given epoch 
reproduces the typical temperature of the WHIM gas in their 
simulations. 

In order to compute the distribution of shocks, we make 
an analogy with the us ual Press-Schechter mass function 
JPress & Schechter 1974). The fundamental assumptions be- 
hind this mass function are: (i) the initial density field is Gaus- 
sian, and (ii) the threshold overdensity for a particle to be 
attached to a halo is given by the dynamics of a collapsing 
spherical top-hat perturbation. In this model, the fraction of 
mass in collapsed halos of mass M > Mq is 

S c (z) 



F(M>M [h z) = edc 



V2a(M ) 



(3) 



Here er(Mo) is the variance in the density power spectrum 
smoothed over a mass scale Mq and S c (z) = 6 c o/D(z) is the 
linearly extrapolated overdensity corresponding to the col- 
lapse of a spherical perturbation. D(z) is the cosmological 
growth factor and <5 c0 « 1 .69 depends very weakly on the as- 
sumed cosmology. The mass function can be obtained by sim- 
ply differentiating equation 0. This mass function matches 
the distribution measured in cosmological simulations to 
good accuracy, although modifications including ellipsoidal 
(rather than spherical) collapse do improve the fit somewhat 
(Sheth & Tormen 1999; Jenkins et al. 2001), though perhaps 
not at the highest redshifts ( Jang-Condell & Hernquist 2001). 

Here we make the ansatz that we can identify shocks with 
a different density threshold S x during the collapse. The en- 
tire Press-Schechter formalism carries over to our situation 
so long as we change the density threshold from that for col- 
lapse to the new value and choose the smoothing scale M s h 
for a specified shock temperature T s h through equations Q 
and 0. The resulting temperature distribution is 
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where df is the fraction of IGM gas (by mass) that is shocked 
to a temperature in the range (r s h, T^ + dT^)- For reference, 
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where h(z) = [O m (l +z) 3 + SIa] 1 ' 2 - Because the T sh -M sh rela- 
tion is itself redshift-dependent, the redshift dependence of a 
cannot be factored out (unlike for the mass function). Also, 
we are interested only in the fraction of gas at a given temper- 
ature, not in the number density of such objects (which would 
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be more sensitive to the non-spherical nature of collapse). In 
all other respects equation (|4j is analogous to mass function 
calculations. 

We note that both 77 and 5 X are essentially free parame- 
ters in our formalism: in the purely spherical model, shocks 
do not actually appear until complete collapse, so the den- 
sity threshold and shock velocity are ad h oc ad ditions to the 
model. However, we will show below (see SI2.1> that a simple, 
physically motivated set of choices provides a reasonable de- 
scription of the existing simulation data for strongly shocked 
material. A plausible point for the shock to first appear is 
turnaround [S x w 1.06/D(z)], the point at which the perturba- 
tion breaks off from the cosmological expansion. In purely 
spherical collapse, of course, no shocks would appear until 
shell-crossing occurs at virialization. However, turnaround 
marks the point at which the flow begins to converge. It is 
therefore reasonable to expect asymmetries in the collapse to 
begin to shock the gas at this point. Another advantage of 
associating shocks with turnaround is that the shock velocity 
in equation ([0 exactly equals the peculiar velocity of the col- 
lapsing gas. Because most shocks are strong, we then assume 
77 = 5/16 for simplicity. We expect this choice to work best for 
regions with infall velocities much larger than the IGM sound 
speed. If the two are comparable, asymmetries in the collapse 
would have to be extremely strong for shocks to develop be- 
fore shell-crossing. In this regime, we might therefore expect 
S x to be closer to S c . 

We also note that our formalism does not describe (or re- 
place) virial shocks around collapsing halos. The "large-scale 
structure shocks" that we describe have much lower velocities 
because they occur earlier in the infall process. If we compare 
the shock temperature r s h on a given mass scale to the virial 
temperature 7^ on the same scale, we find 
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h\z) r) 
so the virial shocks will still be relatively strong. 

There is one more subtlety in computing the temperature 
distribution: a particular trajectory that reaches S x (z) at some 
mass M s h may reach 8 c (z) at another mass Mf, < M sn . Phys- 
ically, the particle could belong to a galaxy embedded inside 
a large-scale filament or sheet. We wish to exclude such gas 
from our distribution because it does not truly belong to the 
IGM. (It will also be at a different temperature: the gas ini- 
tially virializes and then cools.) Fortunately, it is straightfor- 
ward to identify these particles in the Press-Schechter model. 
Given that a particle has (<5 A -,M s h), we wish to compute the 
probability that it has (8 c ,Mh). The problem is formall y iden- 
tical t o the "progenitor" problem solved by Lacev & Cole 
i 19931) . with the solution 

2 a(M h )[S c (z)-S x (z)] 



dp 

~dM h ~ V 7T [a 2 (M h )-a 2 (M sh ,z)]y 2 
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Here dp/dM/, is the probability that a particle is contained in 
a halo of mass M/, at z, given that it is part of a shock with 
mass scale M s h at z- Then 
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We note that this prescription fits nicely into the "merger tree" 
algorithm commonly used in "extended Press-Schechter" 



studies (e.g.. lLacev & Coielll993l) . To compute the probabil- 
ity that a halo of mass Mh will merge into a halo of mass 
M s h at a later time z v , we simply substitute S x (z) — ► S c {z v ) in 
equation |7}. Clearly, z v is fixed by the dynamics of collapse, 
D(z v ) = (S x /5 c )D(z). But this is the same relation as between 
the collapse and turnaround times in the spherical collapse 
model. We therefore find that the probability for a halo to lie 
in a turnaround shock of a region with mass M sn at z is the 
same as for the halo to lie in that region once it has collapsed. 

The correction for bound objects would not be needed if the 
large-scale shocks were able to strip gas from collapsed ob- 
jects. In the local universe, this occurs through ram pressure 
stripping and is fairly inefficient: in galaxy clusters, observa- 
tions and simulation s show tha t the timescale is short only in 
dense cluster cores (Oui lis et alJl20 00). and it is likely to be 
much longer in the low-density IGM. However, if the gas in a 
bound object is unable to cool radiatively and is then shocked 
during infall, it may be able to escape its host halo into the 
IGM. We w ill examine the implications of this process in de- 
tail in JB] 

Note that the postshock temperature depends on the 
whether the gas is ionized or neutral. Any shock with r s h > 
10 4 K will ionize gas after the shock, decreasing the actual 
temperature of the postshock gas by a factor ~ 2. If the gas 
is originally neutral, the shock will also lose energy ioniz- 
ing the gas. For simplicity, we assume that the universe is 
completely ionized for all redshifts lower than the reioniza- 
tion redshift, z < Zr, and neutral for z > z r - The temperatures 
we quote for shocks before reionization do not include ion- 
ization by the shock; they should thus be interpreted as the 
mean thermal energy per baryon. The reionization redshift 
also affects the minimum halo mass able to colllapse. If the 
universe is neutral and cold, we take the "filter mass" Mgi 
or time-ave raged Jeans mass, as the minimum collapse mass 
iGnedin & Huill 998). (Because turnaround shocks occur ear- 
lier in collapse than virialization does, we actually evaluate 
Mm at the redshift corresponding to collapse of the object.) If 
the universe is ionized, we take the minimum mass to have 
r v ; r > 2 x 10 5 K, belo w which gas infall is suppressed by 
photoionization heating ( Efstathiou 1992; Thoul & Weinberg 
1995; Navarro & Steinmetz 1997), although near the time of 
reionizatio n this pro bably overestimates the minimum mass 
( Diikst ra et all2003l) . 

2.1. Comparison to Simulations 

While we have presented a plausible picture in which 77 and 
S x are chosen to match strong shocks at halo turnaround, these 
are actually free parameters that are not fixed by our model. 
In Figure we show how their choice changes our differen- 
tial temperature distribution function, df/dT. The top panel 
shows 77 = 0.16,0.24, and 5/16 (dotted, short-dashed, and 
solid curves, respec tively), with 5 r = 1.06. [Note that the 
simulations of Kesh et etafl (|2003) suggest 77 w 0.16.] The 
sharp rise at low temperatures is the effect of the finite Jeans 
mass that we have imposed. Structures collapsing near the 
mass threshold cannot have any gaseous subcomponents, so 
the fraction of gas associated with bound objects rapidly ap- 
proaches zero. The effect of varying 77 is simply to trans- 
late the distribution along the temperature axis: we have not 
changed the fraction of mass associated with each mass scale, 
only the temperature associated with that scale. The bottom 
panel fixes 77 = 5/16 but sets 6 X = 0.8, 1.06, and 1.20 (long- 
dashed, solid, and dash-dotted curves, respectively). The ef- 
fect here is more subtle: for example, by increasing the thresh- 
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FIG. 1 . — The distribution of gas temperature at z = 0, not including col- 
lapsed objects. Top panel: The curves have r] = 0. 16, 0.24, and 5/16 (dotted, 
short-dashed, and solid lines, respectively); all assume <5 V = 1.06. Bottom 
panel: The curves have & x = 0.8, 1.06, and 1.2 (long-dashed, solid, and dot- 
dashed lines, respectively); all assume r; = 5/16. 



old overdensity required for shocks, we decrease the mass 
fraction at each temperature. Because a(M) is not linear, the 
characteristic temperature also changes (i.e., the peak of the 
temperature distribution). 

Comparison to Figure 5 in iDave et all J2001l> shows that 
our fiducial choice gives a reasonably good fit to the shape of 
the temperature distribution at z = 0. We find that ~ 13% of 
the gas is shocked to high temperature, compared to ^25% 
in the simulations . The comparison is somewhat difficult be- 
cause IDave et all d2001l) include both bound objects (in par- 
ticular galaxy groups) and gas at or near the mean density that 
is nevertheless shocked to high temperatures (see their Fig- 
ure 4), while we expect our formalism to describe only the 
moderately overdense gas. We also find that our formalism 
overestimates the fraction of gas with T s h ~ 10 5 K. This is not 
surprising, because these shocks correspond to systems fairly 
close to the minimum mass for galaxy formation in an ionized 
medium. The infall velocities are therefore not much larger 
than the sound speed of the gas, and shocks probably develop 
only later in the collapse process for such objects. We expect 
our model to apply best when the characteristic shock speed 
is much larger than the sound speed of the IGM. On the other 
hand, capturing weak shocks is difficult in simulations with fi- 
nite resolution dKeshet et alJl2003l iRvu et alJl2003l) . and they 
may underestimate the fraction of gas shocked to moderate 
temperatures. 

2.2. Temperature Evolution of the IGM 

We can now compute the temperature evolution of the IGM, 
including shocks at both turnaround and virialization. The re- 
sults are shown in Figure [2] We assume that the universe is 
reionized at z r = 8; at this point the minimum mass for col- 
lapse increases substantially. The top panel shows the fraction 
of gas that is part of the shocked phase, F^, at a given redshift 
(solid curve) as well as the fraction that is part of collapsed 
objects (dashed curve). The upper and lower dotted curves 
show F s h before reionization if we only include shocks with 
r s h > 10 4 K and 10 5 K, respectively. The bottom panel shows 
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FIG. 2. — Top panel: The fraction of gas in the shocked phase (solid curve) 
and in collapsed objects (long-dashed curve). The upper and lower dotted 
curves show the fraction with T^ > 10 4 K and 10 5 K, respectively. We have 
assumed that reionization occurs instantaneously at z r = 8, at which point 
the minimum halo mass increases dramatically. Bottom panel: The mass- 
averaged temperature of the IGM, including shocked gas only (solid curve) 
and collapsed objects only (dashed curve). The dotted curve shows the tem- 
perature including only shocks with T s h > 10 4 . 



the mass-averaged temperature of the IGM. The solid line is 
for shocked gas while the dashed line is for collapsed gas (the 
temperature of each halo is set to r v i r , and we have used the 
Press-Schechter mass function for the number density of ha- 
los at each mass). The dotted line shows the mean temperature 
from shocks if we include only those with T s h > 10 K (i.e., 
ionizing shocks). We set the temperature of the rest of the gas 
to zero in each case. The thermal energy of the IGM is dom- 
inated by shocked gas at z > 5. At the present time, however, 
virialized objects determine the mean temperature: so long 
as a fair fraction of the mass is in such objects, their strong 
bounding shocks dominate over the IGM shocks [compare to 
equation ©]. 

We see that the fraction in shocks is much larger than that 
in collapsed objects at z > 3, except in a short window near to 
z r - Perhaps surprisingly, the fraction of intergalactic shocked 
gas actually decreases at low redshifts in our model as more 
and more mass is incorporated into collapsed objects. The 
bottom panel of Figure [3] shows the temperature distribution 
function at several redshifts. The solid, long-dashed, short- 
dashed, and dotted lines show results for z = 0, 1, 2, and 3, 
respectively. The top panel shows the fraction of gas nomi- 
nally at each shock temperature that is actually part of bound 
subhalos. The bound fraction rises rapidly above the mini- 
mum mass and then approaches a constant value at a given 
redshift. The asymptotic value is less than unity because of 
the Jeans mass: a non-negligible fraction of the accreted dark 
matter is bound up in potential wells too small to retain gas. 

It is clear that the characteristic temperature of the shocked 
gas declines rapidly with redshift; this is simply because the 
characteristic mass scale also declines rapidly. At z ~ 2-3, 
the nonlinear mass scale is still sufficiently close to the min- 
imum collapse mass that a significant fraction of the gas is 
contained in shocks rather than collapsed objects. By z = 0, 
most of the shocked gas has already entered collapsed ob- 
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FIG. 3. — The distribution of gas temperature in the shocked phase (bottom 
panel) and the fraction of nominally shocked gas that is bound to collapsed 
halos (top panel). Solid, long-dashed, short-dashed, and dotted curves are for 
Z = 0, 1,2, and 3, respectively. All assume S x = 1.06 and r\ — 5/ 16. 
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FIG. 4. — The distribution of gas temperatures in the shocked phase (bottom 
panel) and the fraction of nominally shocked gas that is bound to collapsed 
halos (top panel). Solid, dashed, and dotted curves are for z = 10, 20, and 30, 
respectively. All assume 8 X = 1.06 and 77 = 5/16. 



jects. An important consequence of this is that, at z > 2, 
most of the shocks have T ~ 10 5 K. Because this is close to 
7igm> the typical Mach number of shocks is only ~ 3. We ex- 
pect our prescription for shocks occurring at turnaround to be 
too optimistic in this case, because large asymmetries would 
be required to produce shocks in such a low Mach number 
flow. Thus, our model possibly overestimates the importance 
of s hocks between re ionization and z ~ 2 (compare to Figure 
5 of iDave etal 12001 . Figure[3]demonstrates that (with fixed 
S x and 77) our model does not accurately characterize the im- 
portance of shocks with Mach numbers smaller than 10. The 
fit may be improved by increasing the density threshold at low 
Mach numbers or by decreasing r\ (and hence the shock ve- 
locity). We are primarily interested in the strong shocks prior 
to reionization where this problem does not arise, and so we 
keep the approach simple and fix these two parameters to their 
fiducial values. It should also be noted that the finite resolu- 
tion of the simulations may also cause them to miss shocks (or 
smoot h them out throug h numerical viscosity) at these red- 
shifts; iRvu et ail (12003) find that shocks remain significant 
out to at least z ~ 2, in contrast to the results of Dave et al 
(2001). A more careful treatment of shocks in the simulations 
may provide better agreement with our model. 

We wish to apply our model to the high-redshift IGM in or- 
der to consider the importance of shocks during that epoch. 
The mean temperature of the diffuse IGM is quite small be- 
fore reionization (at least if we neglect X-ray heating). Af- 
ter recombination, the IGM temperature remains coupled to 
the CMB temperature through Thomson scattering off resid- 
ual free electrons. This mechanism becomes inefficient at a 
redshift Zdc given by (Peebles 1993) 

i+z dc «i32(ryi7o.o2) 2 / 5 , (9) 

after which the temperature declines adiabatically as the 
universe expands. This means that 7igm *C 7cmb at z ~ 
20, at least until the first generations of sourc es produce 
suffici ent X-rays to heat the IGM (see 33.21 below and 
Chen & Miralda-Escude 2003). The characteristic shock tem- 
perature, on the other hand, is much larger than this at z < 20. 



Thus infall velocities are large and we expect our model to 
provide a reasonable description. Figure |2 shows that, be- 
fore reionization, the thermal energy of the IGM is dominated 
by large-scale structure shocks rather than virialized objects. 
The disparity gets larger at high redshifts. Essentially, the 
threshold for halo collapse, Mm, is near or even well above 
the nonlinear mass scale at these epochs. The amount of mass 
in collapsing gas is large compared to that in collapsed halos, 
so it dominates the temperature of the IGM. 

Figure |4] shows the temperature distribution for z = 10, 
20, and 30. Again, the characteristic temperature increases 
rapidly with cosmic time. At the highest redshifts, df/dT 
decreases smoothly with temperature and there is no charac- 
teristic temperature peak. This occurs because the nonlinear 
mass scale is below the filter mass at these epochs. This is 
also reflected in the bound fraction. This is small at high red- 
shifts simply because the fraction of gas bound to any halo is 
small. 

Note that these figures do not include radiative cooling. 
While cooling is quite inefficient in the local IGM, it does 
become important at higher redshifts because the characteris- 
tic density increases and the cooling time decreases. Before 
reionization, cooling of regions with T > 10 4 K is quite rapid; 
however, cooling in weaker shocks is negligible because hy- 
drogen remains neutral and the gas is metal-free. Moreover, 
the infalling gas is continually accelerated by gravity, and our 
picture of a single shock at turnaround is most likely too sim- 
ple. In reality, a complex distribution of shocks around and 
inside filaments and sheets may repeatedly shock the infalling 
gas, "refres hing" the energy that has been lost to radiative 
cooling (c.f. lRvu et all200 3 at lower redshifts). We therefore 
expect the curves in the bottom panel of Figure|2]to decrease 
somewhat if radiative cooling is included. In any case, most 
of the applications that we consider are concerned with the 
total energy processed through the shocked gas, which does 
not depend on subsequent cooling mechanisms. 
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3. EFFECTS OF HIGH-REDSHIFT SHOCKS 

3.1. Shock-Accele rated Particles at High-Redshifts 

We first consider cosmic ray acceleration in high- 
redshift shocks. Particle acceleration is know n to be ef- 
ficient in h igh-tempe rature collisionless shocks (Bell 1 1978t 
Blandford & Ostriker 1978). While rare at high redshifts, 
such shocks still contain a substantial amount of energy (see 
Figure^. Even if only a small fraction of this energy is trans- 
ferred to cosmic rays, their energy density can still exceed 
the thermal energy density of the unshocked, neutral IGM. 
We first show that ionizi ng shocks are collisionless, following 
IWaxman & Loebl d2001l) . We will only consider such shocks 
in this section because acceleration is unlikely to be efficient 
if most of the thermal energy is carried by neutral particles. 
For ionized shocks, we wish to compare the ion-ion collision 
frequency, 
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to the growth time for plasma instabilities, v g = C(v s /c)v p i, 
where the plasma frequency is 
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Here C is a factor of order unity that depends on the detailed 
plasma physics. The shock is collisionless if v g > v^, or if 
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We therefore expect par ticle acceleration to be efficient 
( Blandford & Eichler 1987), though we note that the condi- 
tions are sufficiently different from the best studied cases (su- 
pernova blastwaves and cluster shocks) that the properties of 
the acceleration mechanism are not well-constrained observa- 
tionally. 

We assume that plasma instabilities induce a magnetic field 
containing a fraction £g of the thermal energy in equiparti- 
tion. In such a situation, particles are accelerated to hig h en- 
ergies by scattering off fluctuations in the magnetic field ( Bell 
H978t iBlandford & OstrikerlH978l) . The acceleration time is 
face ~ ic/v 2 , where ri is the Larmor radius of the parti- 
cle. The maximum Lorentz factor -f max is fixed by the con- 
dition f acc ~ tic, where fic is the inverse Compton (IC) cool- 
ing time of the particle (the dominant cooling mechan ism at 
high redshifts). T he requirement is ( Loeb & Waxman1l2000t 
Kesh et et all20 03) 
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In fact, even those electrons well below this threshold will 
lose all their energy to IC emission within a Hubble time so 
long as 
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We can therefore assume that all of the energy entering the ac- 
celerated electron component is transformed into IC radiation 



(which is a much more efficient cooling channel than syn- 
chrotron losses at high redshifts given that £b < 1). Assum- 
ing a strong shock, the re sulting spectrum is flat (/„ oc v~ l ; 
Rvbicki & Liahtman 1979). We estimate the energy density 
as in lLoeb & Waxmanl (12000): 
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where T s h mass-averaged temperature of the IGM and £ e 
is the fraction of thermal energy behind the shock that is 
used to accelerate the electrons. The factor 1/2 enters be- 
cause each decade in electron Lorentz factor corresponds to 
two decades in photon frequency. Note that f s h should in- 
clude only those shocks that accelerate electrons (i.e., those 
with r s h > 10 4 K). We show this quantity as the dotted line 
in the bottom panel of Figure |2] Note that, by z ~ 20, 
most of the thermal energy is contained in such shocks. At 
the emission redshift, the IC spectrum will extend between 
(ecMB,7max £ CMB) « (10" 3 , 10 3 ) x (1 +z) eV. The background 
flux is then J v = c/Att x u v : 
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where J v i\ = J„/(10~ 21 erg s" 1 cm" 2 Hz" 1 sr" 1 ) and v a is the 
frequency of the hydrogen Lya transition. 

We thus find that particle acceleration at large-scale struc- 
ture shocks creates a diffuse radiation background with a flat 
spectrum ranging from radio to soft X-ray energies. Interest- 
ingly, such particles may be the first source of high energy 
(> 1 eV) radiation in the universe after cosmological recom- 
bination. We expect particle acceleration to begin once the 
first structures with T s h > 10 4 K turn around. This occurs well 
before z ~ 30 (albeit only rarely). In contrast, simulations 
show that the first stars form at z ~ 20-30 ( lAbel et alJl2002t 
iBromm eta l. 2002); this is both because the star-forming re- 
gions must have fully collapsed and because subsequent cool- 
ing is regulated by H2 formation and is quite slow. We con- 
sider some of the consequences of this radiation background 
in the next section. 

In addition to electrons, shocks will deposit a fraction £ p > 
£ c of their energy in cosmic ray protons. Such protons do 
not cool through IC radiation; in fact, they become an impor- 
tant pressure component of the shock. This component may 
actually have a higher energy density than the diffuse IGM 
(averaged over the entire universe): 
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Note that the shocks continue to accelerate protons until the 
proton Larmor radii exceed the shock thickness. We therefore 
expect the high energy tail of the distribution (which contains 
a substantial fraction of the total energy for a strong shock, see 
equation 11151 ) to escape into the diffuse IGM and fill out the 
voids uniformly. However, the short diffusion time of the rel- 
ativistic ions in the unmagnetized IGM implies that their pres- 
sure has a weak effect on the assembly of gas during galaxy 
formation. 

3.2. Diffuse Radiation 
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One way for the thermal shock energy to be transformed 
into photons is through IC emission by accelerated electrons, 
as described in the previous section. Another is direct radia- 
tive cooling. Figure |2 shows that the thermal energy budget 
of the IGM is dominated by shocked regions before reioniza- 
tion. We therefore expect the background radiation field from 
cooling gas to be dominated by shocked regions as well (pro- 
vided that the cooling time is short in these regions, which 
is true so long as T s h > 10 4 K). Cooling occurs principally 
through IC emission, line transitions, and free-free emission. 
For the moderately overdense regions with T ~ 10 4 — 10 5 K 
in which we are most interested, line cooling typically domi- 
nates. IC emission from thermal electrons produces extremely 
low-energy photons, so we will not consider it further here. 

Most of the cooling occurs through hydrogen and helium 
line emission, which the background spectrum at E < 1 3 .6 eV 
. The flux at a given frequency v is the sum of fluxes from 
those lines that redshift to v at the appropriate cosmic time, 
after adjusting for radiative transfer in the IGM ( Haiman et al. 
(19971) . We can roughly estimate the background flux by 
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where /i; ne is the fraction of the thermal energy lost to this line 
and again r s h should only include ionizing shocks. Note that 
f s h should actually be evaluated at the z' > z appropriate to the 
specified frequency. For example, to evaluate the contribution 
to J v (v a ) at a redshift z from Ly/3 emission, we would set 
\+z' = (v )/(v a )x(l+z). 

Next we can estimate the spectrum of free-free radiation as 
follows. First we approximate the bremsstrahlung emissiv- 
ity n 2 e e(f as flat below a cutoff freqency given by hv c ~ IcT. 
Thus the energy density at a particular frequency is simply 
proportional to the amount of gas that has been shock-heated 
above the corresponding temperature. Each parcel of shocked 
gas radiates over its cooling time (where we must include all 
emission mechanisms). We find 
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where e(T) is the cooling function (including all processes) 
and / s h(> T) is the fraction of gas that has been shock-heated 
to kT > hv. In the first line, the second factor is the ther- 
mal energy in the shocks warm enough to to emit at the rel- 
evant frequency, while the third factor is the fraction of en- 
ergy lost to free-free emission at the appropriate frequency. 
(We have assumed here that the cooling time is smaller than 
the Hubble time; if not, the flux decreases by the appropri- 
ate factor.) The universe is essentially transparent to pho- 
tons with v < v a (corresponding to a temperature threshold 
T S h > 10 5 K), while above this threshold the radiation inter- 
acts with the IGM gas. We show / s h(> 10 5 K) as the lower 
dotted curve in Figure |2 such shocks are extremely rare until 
z ~ 10. We see that IC emission from accelerated particles 
will dominate the high-energy background radiation field at 



X-ray energies, but the free-free emission may be significant 
at UV energies for z < 15. 

While these estimates are quite rough, it is clear that (ne- 
glecting emission from protogalaxies) the radiation spectrum 
from a few eV to the ionization threshold of hydrogen will 
be dominated by free-free line emission from cooling objects, 
with a higher energy tail due to IC emission from accelerated 
electrons. This diffuse background may have dynamical ef- 
fects on structure formation. We list several of these below. 

• Ionizing background: We can estimate the ionization 
fraction due to the IC background as x e ~ /jMic/CifrE,), 
where f\ is the fraction of the radiation background that 
goes into ionizations, is the mean energy lost per 
ionization, and u\c = £ e x (3/2)n b kf s h is the total en- 
ergy density in the radiation field. Note that /} < 1 /2 
accounts for both those photons below the ionization 
threshold and the fraction of energy in high-energy pho- 
tons that goes into heating rather than ionization. The 
resulting electron fraction is 
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The residual electron fraction from cosmological re- 
combination is x e ~ 2 x 10~ 4 with our assumed cosmo- 
logical parameters. We thus find that ionization from 
IC emission could only become important at z < 15 and 
even then only if star formation is quite inefficient and 
hence subdominant. On the other hand, a significant 
fraction of the IC photons have keV energies, and these 
are able to propagate much farther than lower-energy 
stellar ionizing photons. 

• X-ray heating: An additional effect of the high energy 
photons is to heat the diffuse IGM gas. We can estimate 
the resulting temperature by Tx ~ fh£eT S h, where fh is 
the fraction of IC energy that goes into heating the gas. 
Comparing to the temp eratu re of the gas 7igm (without 
additional heating, see ^2.2l above). we find 
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Note that Tx ~ 7cmb at z = 10 with the parameter 
choices shown. Even though only a small fraction of the 
shock energy goes into particle acceleration, the diffuse 
IGM gas is so cold that the heat input can be significant, 
even at z ~ 20. The temperature of the IGM ultimately 
determines the minimum mass of halos that can form 
and is potentially measur able with the 21 cm transition 
of neutral hydrogen (see ^4.2l below). 

• Effects on H2 formation: Without metals, cooling 
to the low temperatu res relevant for s tar formation 



is co ntrolled by H? llAbel et alJ 120021: iBromm et alJ 
2002; Barkana & LoeblhoOll) . Radiative feedback can 
be significant in the formation of molecular hydro- 
gen; ultraviolet photons in the Lyman- Werner band 
(11 .26-13.6 eV) can dissociate H 2 l lHaiman et alJl997l 
iMachacek et~ai1l200l . while X-rays may catalyze its 
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formation by increa sing the density of free electrons 
(Haiman et al. 2000 ), although this effect remains con - 
troversial i lMachacek et afll2003t lOh & HaimarT2 003). 
Photodissociation of H2 becomes important when J21 ^ 
10~ 3 in the UV band; we expect that this could be- 
come significant at z ~ 20 in our model [see equation 
([18}. The IC component has J 2 \ ~ 10" 6 (keV /£), suf- 
ficiently far belo w the UV flux that photodissociation 
should dominate ( Haima n" et all20 00). 

• Wouthuysen-Field Effect: There has been much inter- 
est recently in observing neutral gas at high-reds hifts 
through the 21 cm spin-flip transition of H I (see 34.21 
below). For this to be possible, the spin tempera- 
ture Ts of the hyperfine transition must decouple from 
Tcmb- In the diffuse IGM, d ecoupling occurs th rough 
the Wouthuysen-Field effect (Wouthuvsen 1952; Field 
1958). In this process a hydrogen atom absorbs and 
re-emits a Lya photon; the selection rules allow mix- 
ing of the two hyperfine states during this process. 
The Wouthuysen-Field effect becomes significant when 
A-nOa) > 1 (Ma dauet al.lll997l) . This requires f s 



sh 



> 



10 4 K in equation (II 81 . which does not occur until 
z < 8. We therefore find that, in our model, radia- 
tion from shocks cannot decouple 7^ and 7cmb- Stel- 
lar radiation is required in order to observe the diffuse 
gas through the 21 cm transition (se e ICiardi & M adau 
l2003llChen fc Miralda-Fscudel2003h. 

Of course, both IC and cooling radiation come from dis- 
crete sources and are not truly diffuse backgrounds. This is es- 
pecially true at high redshifts, where shocks with T s h > 10 4 K 
are quite rare. We therefore expect the above effects to be 
significantly larger near to the shocks themselves. Ioniza- 
tion and feedback effects on H2 would be considerably larger 
in these regions, which is especially important because both 
shocks and collapsing structures preferentially form in the 
most highly biased regions. In principle, the cooling radiation 
from individual shocks could be imaged. Unfortunately, while 
the luminosities are large, the sources are quite extended and 
the surface brightness is extremely small. Even the James 
Webb Space Telescope would require hundreds of hours of in- 
tegration to see an individual system at z > 10. 

3.3. Gas Stripping from Minihalos 

At low and moderate redshi fts, large-scale s tructure shocks 
have complex morphologies (iRvu et alJl2003l) that surround 
and interpenetrate sheets and filaments. Nearly all of the gas 
that accretes smoothly onto the structures is shocked; how- 
ever, it is less clear whether gas that has cooled onto galaxies 
or even galaxy groups is affected by the shocks. The galac- 
tic interstellar medium is much denser than the surrounding 
IGM, which dramatically decreases the strength of the accre- 
tion shocks felt by the interstellar medium. Moreover, these 
objects have substantial magnetic fields that may deflect gas 
in the external medium around the object, rather than allowing 
it to collide and shock the interstellar medium of the galaxy. 

However, these mechanisms do not apply to gas in "mini- 
halos" at high redshifts. Such objects have virial tempera- 
tures r v i, < 10 4 K, placing them below the threshold to ex- 
cite atomic hydrogen transitions. In the absence of molec- 
ular hydrogen, they are therefore unable to cool efficiently 
and cannot collapse after virialization, remaining as moder- 
ately overdense (S ~ 200) gaseous halos in the IGM (see 



iBarkana & Loeblf2001l and references therein). In terms of 
their internal structure, such halos are thus more similar to 
local galaxy clusters than to galaxies, except that they likely 
have negligible magnetic fields and no gas substructure. 

Thus, these halos have little protection from the ram pres- 
sure of the surrounding medium. The minihalo gas will 
then be stripped to a radius R determined by the condition 
p sh v* h = PmhiR)^, where p sh = S sh p is the mean density in 
the shocked region, p m h(R) = S m h(R)p is the density at a given 
point within the minihalo, and cr m h is the velocity dispersion 
of the minihalo. The condition is approximately 
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Given a density profile for the minihalo gas, this prescription 
allows us to compute the fraction of minihalo gas that is lost 
as it travels through large-scale structure shocks. We consider 
two possible density profiles: a singular isothermal sphere 
(SIS) and isothermal gas inside of a dark matter halo with 
the structure found by Navarr oet al.l 111997. hereafter NFW). 
We parameterize our results in terms of / ret , the fraction of the 
initial minihalo gas mass retained by the minihalo after ram 
pressure is taken into account. For an SIS, the criterion is par- 
ticularly simple: / let = R/r wlI . For isothermal gas in an NFW 
halo, 

/ret W = J °„ 1 ; ; — — , (23) 
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and c is the concentration parameter of the minihal o. We 
choose c from the fit to simulations given by B ullock et alJ 
(2001). 

Figure |3 shows the effects on individual halos. The bot- 
tom panel shows the probability that a halo at z = 10 is con- 
tained within a shock as a function of T^/T^. The solid, 
long-dashed, short-dashed, and dotted curves assume r v ; r = 
170, 1080, 5030, and 10 5 K, respectively. These correspond 
to M h = 6 x 10 4 , 10 6 , 10 7 , and 10 9 M Q . The smallest of the 
these masses is the filter mass at z = 10 while the largest is 
above the cooling threshold and is included only for illustra- 
tive purposes. We see that the probability for a halo to sit 
inside a strong shock is extremely sensitive to its mass. A 
substantial fraction of halos below the nonlinear mass scale 
are embedded in strong shocks, but very few halos above that 
scale are associated with larger shocks because such struc- 
tures are increasingly rare. The top panel shows the fraction 
of mass retained in a halo as a function of T^/T w - U . The solid 
line assumes the gas has an SIS distribution (/ ret is then inde- 
pendent of halo mass when expressed as a function of 7^/7^) 
and the dot-dashed lines assume that the gas is isothermal in 
an NFW halo. The concentration parameters in this plot are 
for r v j, = 170 and 10 5 K (left and right curves); this is the only 
aspect in which the halo mass enters into this parameteriza- 
tion. Stripping sets in when the shock temperature is a few 
times the virial temperature of the minihalo. The SIS model 
is more concentrated and is therefore more stable to mass loss. 
Note that the NFW model has a finite density core, so the en- 
tire minihalo can be stripped. Comparing the two panels, we 
see that a negligible fraction of mass would be lost from halos 
near to the minihalo mass threshold, but a large fraction could 
be lost from smaller halos. 
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FIG. 5. — Bottom panel: The probability that a halo at z = 10 is contained 
within a shock of a given strength. The solid, long-dashed, short-dashed, and 
dotted lines are for halos with T vir = 170, 1080, 5030, and 10 5 K. Top panel: 
The fraction of mass retained in a halo after ram-pressure stripping. The solid 
line assumes the gas has an SIS distribution and the dot-dashed lines assume 
that the gas is isothermal in an NFW halo. In the latter case, the left and right 
curves assume concentration parameters appropriate for halos with r v j, = 170 
and 10 s K, respectively. 



Figure [5] shows the global effects of shock disruption. We 
show the fraction of mass in collapsed objects (including both 
minihalos and halos that are able to cool) that remains after 
ram pressure stripping is included. The solid line assumes 
that minihalos have SIS density profiles while the dashed line 
assumes they are isothermal gas in an NFW potential. The 
two density profiles give very similar results because they dif- 
fer substantially only in the core (which rarely gets stripped). 
We see that the global effects of stripping are fairly small: 
on ly ^ 5% of the gas is stripped at z > 20 while approxi- 
mately 10% of the gas is stripped at z ~ 10. The importance 
of stripping increases with cosmic time because, as the non- 
linear mass scale increases, regions with stronger and stronger 
shocks form and a larger fraction of minihalos are contained 
in such regions. Note that stripping is concentrated in the most 
highly biased regions of the universe (where strong shocks 
can form). This is also where a large fraction of ionizing 
sources reside; we would therefore expect minihalos to be less 
important aroun d such sources than previous estimates hav e 
assumed (e.g., Barkana & Loeb 2002; Shapiro et al. 2003). 
These estimates also assume that minihalo gas is stripped only 
after the minihalo has formed. In reality, many minihalos will 
enter large-scale structure shocks during their own collapse. 
Criterion \22\ suggests that stripping is more efficient if the 
initial density is smaller. 

Ram pressure is not the only way in which shocks can strip 
minihalos. As they fall into these structures, the minihalo gas 
will pass through the shock. If the shock adds suffici ent en- 
tropy to the gas, it can escape the minihalo ( Oh & Haiman 
2003). The stripping efficiency is then determined by the ra- 
tio of the shock speed in the minihalo [which is proportional 
to p(R)~ l l 2 ] to the sound speed. We find that this mecha- 
nism is slighly less efficient than ram pressure stripping, so 
the two mechanisms together will likely increase the amount 
of stripped gas over that shown in Figure [6] by a small fac- 
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FIG. 6. — Fraction of gas in collapsed objects that is retained after ram 
pressure stripping. The solid curve assumes the minihalo profile is an SIS, 
while the dashed curve assumes isothermal gas in an NFW potential. 



tor. Calibration of the precise level of these effects must await 
numerical simulations that include a careful treatment of the 
shock physics as well as broad enough dymanic range to re- 
solve minihalos while also encompassing a large, representa- 
tive cosmological volume. 

Finally, note that shocks may trigger star formation in mini- 
halos. If the postshock temperature > 10 4 K, the gas 
will be able to cool through both line radiation and inverse- 
Compton processes. In this case, gas would remain bound 
to the minihalo and, if cooling is sufficiently fast, some gas 
could even collapse to hi gh enough densities to form stars 
(Mackev et al. 2003; Cen 2003). In reality, the dense central 
regions of the minihalos may collapse to form stars while the 
outer layers are stripped by the shock. Disentangling these 
two effects requires hydrodynamic simulations of the shock- 
ing process. Triggered star formation is likely to be less effi- 
cient in large-scale structure shocks than in supernova blast- 
waves because the characteristic velocities are smaller. 

4. OBSERVATIONAL SIGNATURES OF HIGH-REDSHIFT SHOCKS 

4.1. Optical Depth to the CMB 

One way to observe the effects of gravitational heating 
is through its effects on the CMB. Large-scale shocks with 
T > 10 4 K will ionize gas that passes through them. CMB 
photons will scatter off the resulting free electrons, smoothing 
out the anisotropies. The cumulative optical depth to electron 
scattering r es between redshifts zi and Zh is 

f~" dt 
r es = J dza T n e (z)c—^ (25) 

where aj is the Thomson scattering cross section. We com- 
pute the optical depth by assuming that all gas associated with 
shocks of T > 10 4 Kat a given redshift is ionized. The result 
is shown in the top panel of Figure0 T es < 1 % and is unlikely 
to contribute substantially to the measured optical depth, par- 
ticularly if reionization occurs at high redshift. By neglecting 
recombination of the gas, we actually provide an upper limit 
to the optical depth due to IGM shocks. Recombination will 
be significant within the overdense shocked regions, unless 
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FIG. 7. — Top panel: The optical depth to electron scattering r cs in our 
model, assuming that all gas associated with bound objects and/or turnaround 
shocks of T > 10 4 K remains ionized. Bottom panel: The mean 21 cm bright- 
ness temperature decrement ST[> of the IGM relative to the CMB. The dotted 
curve neglects all heat sources (including gravity). The dot-dashed curve in- 
cludes virial shocks. The solid and long-dashed curves include heating at 
turnaround shocks; the latter assumes that all gas in shocks with T > 10 4 K 
remains ionized and is almost hidden by the solid curve. The short-dashed 
curve also includes X-ray heating through the IC background from shock- 
accelerated electrons. 



the gas is shocked many times as it falls into the potential 
well of an object. 



4.2. 21 cm Radiation 



Another way to identify shocks is through their effect on 
the global temperature distribution of the IGM, which can be 
measured with the 21 cm transition. The brightness tempera- 
ture increment of the 21 cm transition from neu tral gas relative 
to the CMB is ( Fie ld! 1 95 8tlMadau et alJl997i) 
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where T s is the spin temperature of the hydrogen and the inte- 
gral extends over all neutral gas. Note that if 7^ 3> 7cmb, the 
signal saturates in emission, but if Ts <§C 7cmb the magnitude 
of the absorption increases like T$ l . 

The spin temperature remains locked to the CMB temper- 
ature unless (1) the background Lya photon field is large 
enough for t he W outhuysen-Field mechanism to become ef- 
ficient (see 33.21 above) or (2) the number density of hydro- 
gen is large enough for collisions to decouple the two ( Field 
1958). If these processes operate, 
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where T a is the Lya color temp erature (T u = Tk under a wide 
variety of conditions: lFieldl 1 9591) and Tk is the kinetic temper- 
ature of the gas. The coupling constant y a = Pi<jT+/(An)T a ), 
where P w = 1.3 x \0~ l2 J Vj ii(v a ) is the indirect de-excitation 
rate of the triplet level due to the absorption of a Lya pho- 
ton followed by decay to the singlet level, = 0.068 K is 



the excitation temperature of the hyperfine level, and A\o = 
2.85 x 10~ 15 s" 1 is the spontaneous emission coefficient of the 
transition. The coupling from collisions is y c = C\oT+/ (Au)T a ), 
where Cm is the collisional de-excitation rate of the triplet 
hyper fine level and has b een calculated for Tk < 10 3 K by 
Allison & Dalgarno ( 1969) (we assume Cio oc 7^ 33 for higher 
temperatures). 

We first assume that T$ = Tk throughout the unive rse ( i.e., 
there is a strong Lya background). We showed in 33.21 that 
this requires stellar (or quasar) radiation, because the back- 
ground from shocks themselves is insufficient. In the absence 
of heat sources, the IGM will appear as a smooth absorption 
feature against the CMB [df/dTs = S(T -Tigm)]- The signal is 
shown by the dotted curve in the bottom panel of Figure0 we 
see \5Th\ > 100 mK, decreasing rapidly as redshift decreases 
because Tk oc (1 +z) 2 . This signal is sufficiently strong to be 
easily observable by the next generation of low-frequency ra- 
dio telescopes. However, structure formation will remove gas 
from the cold IGM, decreasing the overall absorption signal. 
First, we must remove gas that is attached to bound objects. 
In this case the integral in equation (1261 extends over the mass 
function of collapsed objects dn/dM. The result is shown in 
FigureQas the dot-dashed line. We see that ST/, changes by 
only a small amount for z > 15. 

Next, we include (in addition to collapsed objects) gas that 
has been shock-heated by large-scale structure. This is shown 
by the solid line in Figure assuming that all gas remains 
neutral. (The long-dashed line shows the signal if gas asso- 
ciated with T > 10 4 K shocks remains ionized and does not 
emit; it is nearly completely obscured by the solid curve.) 
Because / s h / co ii at high redshifts, the absorption signal 
flattens out at much higher redshifts; even at z = 20, |<57/,| is 
smaller by about 30 mK. In principle, shock-heating at high 
redshifts could therefore be detected by its effects on the radio 
spectrum. However, other heat sources, such as X-rays from 
the IC back grou nd generated by particles accelerated at the 
shocks (see 43. 2\ or from supernovae or quasars, will increase 
the temperature of the low-density IGM and may cause signif- 
icant contamination with this signal llChen & Miralda- Escudel 
2003). For example, the short-dashed line shows the signal 
assuming the heat input from equation i2l\ . which includes 
IC emission from shock-accelerated electrons. Even this rela- 
tively gentle heating rapidly decreases the amount of absorp- 
tion, and it actually turns into emission by z ~ 13. We would 
therefore require an early stellar Lya background (unaccom- 
panied by X-ray heating) for significant absorption to be pos- 
sible. Note that the treatment of weak shocks has a substantial 
effect on this estima te, so our model may not be quantitatively 
accurate. Gnedin & Shaver (2003) have recently simulated 
this absorption era; they find that shock heating decreases 
\STt,\ by an even larger amount than we estimate. This in- 
dicates that weak shocks preceding turnaround may heat the 
diffuse IGM by a small amount before reionization and that 
we underestimate their effects (unlike at moderate redshifts). 

The above discussion of the globally-averaged signal ig- 
nores the fact that the temperature and density distributions 
are highly inhomogeneous. The absorption and emission sig- 
nal therefore fluctuate s across the sky, and one can also search 
for these fluctuations dScott & ReesH990tlMadau et all 19971 
Zaldarriag a et all2003U Loeb & Zaldarriagga 2003). Because 
the emission signal is independent of Ts for Ts 3> 7cmb, it 
essentially depends only on the mass of gas in a given vol- 
ume that emits. Let us now consider a sufficiently early epoch 
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when the Lya backgrou nd is s mall, so that Ts = 7cmb in the 
diffuse IGM. H ievetal] (12002) have shown that, because col- 
lisional coupling is still effective in dense gas, minihalos pro- 
vide a fluctuating emissi on signal t hat may be observable on 
large scales (see also Iliev et al. 2003). Here we point out that, 
because the shocked phase contains so much more mass than 
minihalos, the e mission signal is probably much larger than 
llliev et al\ (120021) have estimated. 

In the case of minihalos, the density is so large that y c 3> 1 
throughout the halo. Thus T$ = Tg is a good approximation 
in these systems. In shocks, however, the densities are much 
smaller and not all of the gas may have 7$ > Jcmb- Given z 
and Tk, we first define #21, the overdensity for which y c = 1 
and above which collisional coupling becomes efficient. Note 
that (at fixed temperature) 621 increases as z decreases because 
the mean density also decreases. We then let (S21) be the crit- 
ical overdensity for collisional coupling to be effective, aver- 
aged over the temperature distribution of the shocked gas: 

(S 2l (z)) = J dT K ^-5 2 i(z,T K ) (28) 

The resulting values are shown by the solid line in the top 
panel of Figure|S]for our shock model. The dotted line shows 
the overdensity at turnaround in the spherical model; at z > 
25, all shocked gas emits strongly. Note that (<52i) <C 200 for 
all redshifts shown. 

We wish to find the fraction of collapsing gas that has 
5 > In the spherical collapse model, the overdensity 

at turnaround is <5 ta = 5.55 a nd increases m onotonically with 
time until virialization (e.g. Peeblesl H980l) . We can estimate 
the relevant fraction with a slight extension of our model: we 
assume that all gas collapses spherically, so that a parcel of 
gas with a given density at z has a well-defined turnaround 
redshift z ta . The fraction of gas with 5 > (621) at z is then 
simply the fraction of gas that has passed through turnaround 
shocks by Zm- Note that our model is only approximate be- 
cause we do not follow (<5, Ts) of individual gas elements; we 
actually underestimate the fraction of gas with Ts > Tk be- 
cause the hottest shocks, associated with the most massive 
objects, have smaller critical densities and hence may have 
collapsed later on. The bottom panel of Figure [8] compares 
the fraction of the IGM gas that is able to emit 21 cm ra- 
diation in our shock model (solid line) with the fraction of 
gas in minihalos (dashed line). The peak at z ~ 25 occurs 
because (621) < at higher redshifts, so that all shocked 
gas was able to emit. At z > 15, the former is much larger 
than the amount of mass in minihalos. Thus, we see explic- 
itly that the true emission signal will be significantly larger 
than that calculated just from minihalos. Weak shocks occur- 
ring before turnaround ( as appear to exist in the simulations of 
iGnedin & Shaverl200 3) would further amplify the fluctuation 
signal. 3 

5. DISCUSSION 

In this paper, we have presented a model for the temperature 
distribution of IGM gas that takes into account shocks during 
structure formation. We use a simple extension of the Press- 
Schechter model that associates shocks with a fixed density 
contrast, and we assume that the shock temperature is deter- 
mined by the infall velocity at the shock time. Our model 

3 Note that the bias amplification for shocked gas is somewhat smaller than 
that for minihalos, because it is associated with smaller density peaks. How- 
ever, on the large scales accessible to the next generation of radio telescopes, 
the difference should not be large. 
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FIG. 8. — Top panel: The critical density (<52i) for coupling between 
the spin and kinetic temperatures of H I, averaged over the temperature of 
the shocked phase. The dotted line shows 8 = 5.55, the nominal density at 
turnaround. Bottom panel: The fraction of gas able to emit. The solid line is 
for the shocked phase (see text); the dashed line includes only the collapsed 
objects. 



therefore describes moderately overdense gas in the process 
of collapsing. We find that associating shocks with turnaround 
adequately describes the temper ature d istribution of shocked 
gas at z = found bv lDave et alJ OOOl). The agreement with 
simulations is not good for weak shocks (which dominate at 
z ~ 2-3), but we are primarily concerned with strong shocks. 
We then used our formalism to consider shocks at higher red- 
shifts. We found that large-scale structure (rather than viri- 
alization) shocks dominate the thermal energy content of the 
IGM for most of the period before reionization. As a result, 
shocks can have several interesting consequences for struc- 
ture formation. These include stripping gas from ~ 10% of 
all minihalos, depositing a substantial fraction of their energy 
in cosmic rays, and creating a radiation background (due to 
cooling radiation and inverse Compton emission by acceler- 
ated electrons) that extends to high energies. We have argued 
that these shocks can be observed through their effects (both 
directly and indirectly) on redshifted 21 cm fluctuations. The 
warm, moderately overdense shocks will emit 21 cm radia- 
tion and may be visible as a fluctuating background. The 
fluctuations will be larger than those from collapsed objects 
because significantly more gas is contained in the shocked 
phase. Moreover, the X-ray background produced by shock- 
accelerated electrons significantly affects the temperature of 
gas in the low-density IGM, decreasing (or even eliminating) 
any absorption signal. 

Our formalism involves several important simplifications. 
First, we assume that infalling gas is shocked only twice: once 
at turnaround and once at virialization. We have also assumed 
pure spherical coll apse. This is clearly only an approxima- 
tion; iKeshet et alJ ( 120031) and lRvu etatl ( 120031) have shown 
that the distribution of shocks is much more complicated, with 
strong shocks bounding filaments and sheets with a network 
of weaker shocks surrounding virialized objects. The non- 
spherical nature of collapse is also critical, because it affects 
both the fraction of gas that is shocked (through the density 
threshold required before significant collapse occurs) and the 
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strength of the resulting shocks (because the infall velocities 
also change). 

Another caveat is that our model associates shocks with 
turnaround. On the one hand, we seem to overestimate the 
fraction of gas that is weakly shocked at moderate redshifts. 
We have argued that this is because gas flows are only moder- 
ately supersonic in such shocks. When infall velocities are 
slow, it is plausible that shocks can only form later in the 
collapse process. On the other hand, if gas travels super- 
sonically during the early stages of infall, asymmetries in the 
flow may cause parcels of gas to collide and shock well be- 
fore turnaround. This is particularly important before reion- 
ization, when the IGM gas temperature is ~ 10 K. In this 
regime, supersonic flow occurs well before turnaround and 
a much larger fraction of the IGM may be weakly shocked. 
However, the thermal energy of the shocks is dominated by 
the largest objects if z < 20 (compare the solid and dotted 
curves in Figure[2). The strong shocks at or after turnaround 
therefore probably contain most of the energy. Furthermore, 
most of the effects that we consider at high redshifts (includ- 
ing cooling radiation, particle acceleration, and ram pressure 
stripping) depend primarily on the characteristics of the strong 
shocks that are able to ionized the cold pre-reionization IGM. 
We therefore expect our conclusions to be robust against un- 
certainties in the treatment of weak shocks. Weak shocks will, 
however, substantially increase the 21 cm signal during the 
early phase s of structure formation (compare our FigureQand 
Figure 2 of lGnedin & Shave J2003ft . 

Fina lly, we have assumed that the iPress & Sch echter 
J1974I) formalism prov ides a good fit to the m ass function. 
ISheth & Tormenl i ll 9991) and lJenkins et alJ d200 ll) have shown 
that a modified version of the formalism incorporating el- 
lipsoidal collapse provides a better fit to many simulations 



(but see iJang- Condell & Hernauistll200l . We have chosen 
not to use these modifications because we require the condi- 
tional mass functions to remove mass attached to collapsed 
objects. Conditional functions can be written analyti cally i n 
the standard Press-Schechter approach tlLacev & Coldll993h . 
The Sheth-Tormen approach, however, is based on a "mov- 
ing barrier" for which conditional mass functions cannot be 
solved analytically, al though they can be fit to simulation data 
(Sheth & Tormen 2002). For our purposes, the major differ- 
ence between the two approaches is that highly nonlinear ob- 
jects are more common in the Sheth-Tormen mass function. 
We therefore expect high-temperature shocks to be somewhat 
more important than our model predicts. 

Obviously, the best way to test our model is with high- 
resolution simulations including a careful treatment of shocks. 
Such simulations would follow the non-spherical nature of the 
collapse directly and can (at least in principle) capture the 
shocks at many different stages of collapse. However, note 
that the most significant shocks at high redshifts (those that 
ionize the infalling gas) are quite rare. Large-volume simu- 
lations are required in order to include a re presentative sam- 
ple o f these rare and highly biased peaks ( Bar kana & Loebl 
2003). Shock-capturing schemes could then be used in these 
simulations to predict more accurately the background radi- 
ation from accelerated particles, the temperature distribution 
and its effects 21 cm signal, and the importance of ram pres- 
sure stripping. 
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